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^ '. ABSTRACT 

' We combine constraints on the galaxy-dark matter connection with structural and dy- 

, namical scaling relations to investigate the angular momentum content of disc galax- 

ies. For haloes with masses in the interval 1Q 11S M@ < Af vir < 10 12 ' 7 M Q we find 
r*"** ' that the galaxy spin parameters are basically independent of halo mass with (A' ,) = 

, ,' (Jgal/M g ai)/(V^i2virKir) = 0.019±g:gg| (Ict"). This is significantly lower than for re- 

, laxed ACDM haloes, which have an average spin parameter (A^ alo ) = 0.031 ± 0.001. 

r \ • The average ratio between the specific angular momentum of disc galaxies and their 

^■^ ' host dark matter haloes is therefore IZj = A' j/Au alo = 0.61^q \{. This calls into ques- 

■ tion a standard assumption made in the majority of all (semi-analytical) models for 

(disc) galaxy formation, namely that IZj — 1. Using simple disc formation models 

q . we show that it is particularly challenging to understand why IZj is independent of 

!h ' halo mass, while the galaxy formation efficiency (cqf, proportional to the ratio of 

c/3 \ galaxy mass to halo mass) reveals a strong halo mass dependence. We argue that 

C$ . the empirical scaling relations between cqf, TZj and halo mass require both feedback 

(i.e., galactic outflows) and angular momentum transfer from the baryons to the dark 

\ matter (i.e., dynamical friction). Most importantly, the efficiency of angular momen- 

■ turn loss need to decrease with increasing halo mass. Such a mass dependence may 
^C) ' reflect a bias against forming stable discs in high mass, low spin haloes or a transition 

, from cold-mode accretion in low mass haloes to hot-mode accretion at the massive 

• end. However, current hydrodynamical simulations of galaxy formation, which should 

include these processes, seem unable to reproduce the empirical relation between €gf 

qq . and IZj . We conclude that the angular momentum build-up of galactic discs remains 

f^i ■ poorly understood. 

^ | | Key words: galaxies: formation — galaxies: fundamental parameters — galaxies: 

• • , spirals — galaxies: structure 
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1 INTRODUCTION make galaxy discs with the observed sizes (e.g., Dalcanton 

, , „ , r , , et al. 1997; Mo, Mao, & White 1998; Firmani & Avila-Reese 

In the current paradigm of galaxy formation, galaxy discs ™™ it d t „„„„ AT c ni ■ , 

. , . , r i, 2000; de Jong & Lacey 2000; Navarro & Stemmetz 2000). 

are considered to form from the accretion of gas inside hier- T r , , , . , . , , . , o^m- i 

° in fact, models m which galaxies have ~ 60% lower spe- 



archically growing cold dark matter (CDM) haloes (see Mo 
van den Bosch & White 2010 for a comprehensive overview). 
The dark matter (DM) and gas acquire angular momentum 



cific angular momentum than their dark matter haloes may 
be preferred (Dutton et al. 2007). Another success of this 
picture is that the large variation in the specific angular 

cx- 



via tidal torques in the early Universe (Peebles 1969). When c , , , , , , , , 

, , ' momentum of dark matter haloes, at fixed halo mass, 

gas accretes onto the central galaxy, this angular momentum , . . . , , r 

plains the large scatter m disc sizes, or equivalently surface 
may eventually halt the collapse and lead to the formation of , . , , . „ , , . , , , -. 

. „ ... ,„ „ „. , . ,„„„\ T r , brightnesses, at hxed luminosity (Dalcanton et al. 1997). 

a rotationally supported disc (lall & rjjfstatmou 1980). If the ™ . . , , , „ ,. , ,. r ,. . 

, r , . ' 1 his picture also accounts for the time evolution of disc sizes 

specihc angular momentum of the baryons is conserved dur- ,„ .,, , . „„„„ . „ . ., „ nnr , n ^ 

, . , , , (Somerville et al. 2008; Firmani & Avila-Rccsc 2009; Dut- 

mg galaxy formation, there is enough angular momentum to , < , . , , , . ,. 

ton et al. 2011a), reproducing the observed weak evolution 

since redshift z ~ 1 (Barden et al. 2005), and the stronger 
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evolution observed from z ~ 2 to z ~ 1 (Trujillo et al. 2006; 
Williams et al. 2010). 

Despite the many successes of these models, all of which 
assume that specific angular momentum is conserved during 
galaxy formation, there are several reasons, from a physical 
point of view, as to why the ratio between the specific angu- 
lar momentum of the baryons and that of the dark matter 
halo, denoted 7?.jQ, might be different from unity: 

(i) Selective accretion of baryons. If not all available 
baryons make it into the galaxy, the specific angular mo- 
mentum of the galaxy can be either higher (TZj > 1) or lower 
(TZj < 1) than that of all baryons (whether it is higher or 
lower simply depends on the specific angular momentum dis- 
tribution of the materials that are selected) . An example of 
a selective accretion process is cooling. Cooling is generally 
an inside-out process, with the inner regions cooling prior to 
the more outer regions. Since the mass of ACDM haloes is 
more concentrated than the specific angular momentum (see 
e.g., Navarro & Steinmetz 1997), one expects that TZj < 1 
if not all baryons manage to cool. 

(ii) Exchange of angular momentum between the dark 
matter and the baryons. An example of this is dynamical 
friction, which transfers orbital angular momentum from the 
baryons to the dark matter, thus resulting in TZj < 1 (e.g., 
Navarro & Benz 1991; Navarro & White 1994). The strength 
of this process depends on the efficiency of star formation at 
early times. If the formation of baryonic clumps can be sup- 
pressed, such as through strong feedback, then the angular 
momentum transfer is reduced (e.g., Weil, Eke & Efstathiou 
1998; Sommer-Larsen, Gelato & Vedel 1999; Eke, Efstathiou 
& Wright 2000). 

(iii) Feedback. If feedback (from supernova, AGN or any 
other mechanism) results in the removal of gas, one can 
achieve values for TZj that are very different from unity 
by simply removing a specific subset of the baryons. As an 
example, several studies have argued that supernova feed- 
back affects low angular momentum material more than high 
angular momentum material, causing a net increase of TZj 
(Mailer & Dekel 2002; Dutton & van den Bosch 2009; Dut- 
ton 2009; Governato et al. 2010; Brook et al. 2011). 

Cosmological hydrodynamical simulations almost in- 
variably predict that disc galaxies have significantly less 
specific angular momentum than their dark matter haloes 
(TZj < 1), suggesting that mechanisms (i) and (ii) dis- 
cussed above play an important role. Early simulations 
in a standard cold dark matter (SCDM) cosmology (i.e. 
Q m = 1,0a = 0) resulted in a catastrophic loss of angu- 
lar momentum (Navarro & Benz 1991; Navarro & White 
1994; Steinmetz & Navarro 1999), producing discs an order 
of magnitude too small. More recent simulations in A Cold 
Dark Matter (ACDM) cosmologies have resulted in less an- 
gular momentum loss, but typical ratios between the spe- 
cific angular momentum of the baryons and that of the dark 
matter are TZj ~ 50% (Piontek & Steinmetz 2011; Sales 
et al. 2010). This "angular momentum loss" is sometimes 



1 Note that in the notation of Mo, Mao, & White (1998), TZj = 
(id/ m d)i where is the ratio between the angular momentum 
of the disc and dark matter halo, and mj is the ratio between the 
mass of the disc and dark matter halo. 



attributed to lack of numerical resolution (e.g. Governato 
et al. 2004; Kaufmann et al. 2007). However, Piontek & 
Steinmetz (2009) argue that the dominant cause is dynami- 
cal friction between clumps of baryons and the dark matter 
halo. If, based on these findings, we assume that angular 
momentum loss is a natural feature of disc galaxy formation 
in a ACDM universe, a relevant question is the following: 
If the baryons that make up a discy galaxy carry half the 
specific angular momentum of the dark matter halo in which 
they reside, is it possible to reproduce in detail the observed 
sizes of disc galaxies? 

This question was addressed by Navarro & Steinmetz 
(2000). Using scaling relations between disc size, luminos- 
ity, and rotation velocity, these authors derived an expres- 
sion for TZj (fj in their notation) that depends on the galaxy 
formation efficiency, e G F = (M sa ,i/M vil )/(D. h /Q m ) (/bdsk in 
their notation), cosmic matter density, Q m , and stellar mass- 
to-light ratio, Yj. They found that if the rotation speeds of 
galaxy discs are approximately the same as the circular ve- 
locities of their surrounding haloes (which has subsequently 
been shown to be the case, at least on average, by Dutton 
et al. 2010b), then discs must have retained about half of the 
available specific angular momentum during their assembly. 
Another result from this analysis is that disc galaxies form- 
ing in low density cosmogonies, such as ACDM, need to 
have TZj ~ 2 cgf, that is, they draw a larger fraction of the 
specific angular momentum than the galaxy formation effi- 
ciency. This is at odds with the expectation (see § l4.1l below) 
that TZj ~ £gf if the galaxy formation efficiency is mainly 
determined by the inefficiency of cooling. 

Thus observational measurements of the specific an- 
gular momentum in galaxies place constraints on the way 
galaxies acquire their mass and angular momentum. The 
specific angular momentum of the dark matter, J v ir/M v i r , 
where J v i r is the total angular momentum, and M v i r is the 
total virial mass, is commonly expressed in terms of the di- 
mensionless spin parameter (Peebles 1969): 



J vir |£| 1/2 _ J vir /M 



GM 



5/3 



a/2 i^vir Vvi 



vir »l/2 _ w rl/2 
Jc ^halo Jc 



(1) 



Here E is the halo's energy, i? v i r is the virial radius, V v i T is 
the circular velocity at the virial radius, and f c measures 
the deviation of E from that of a singular isothermal sphere 
truncated at i? v ir. It is common to set f c = 1 (Bullock 
et al. 2001), so that the spin parameter only depends on 
the total angular momentum and mass of the halo. In this 
case the spin parameter is denoted A halo . The equivalent spin 
parameter for the galaxy is defined as 



(2) 



(see Section [2. II for details). 

Dark matter only cosmological simulations have shown 
that both Ahaio and A halo are log-normally distributed (e.g., 
Bullock et al. 2001) with median and scatter independent 
of halo mass (e.g., Maccio et al. 2007; Bett et al. 2007; 
Mufioz-Cuartas et al. 2011) and independent of environ- 
ment (Lemson & Kauffmann 1999, Maccio et al. 2007). 
For relaxed haloes identified with a spherical overdensity 
(SO) algorithm Maccio et al. (2008) find a median value of 
Ahaio = 0.031 ± 0.001. Using a similar halo definition Bett 
et al. (2007) find a median Ahaio = 0.0367. For the cosmol- 



© 2012 RAS, MNRAS OOO.rflfTsI 



Angular momentum of disc galaxies 3 



ogy adopted by Bett et al. , f c ~ 1.3 for a 10 M© halo, 
and thus A halo ~ 0.032, which is in excellent agreement with 
the results of Maccio et al. (2008). Simulations with gas in 
which there is no cooling find similar spin parameters for the 
baryons and dark matter (van den Bosch et al. 2002; Sharma 
& Steinmetz 2005). Thus in what follows we assume that the 
median spin parameter of the halo is 

Ahaio = 0.031 ± 0.001. (3) 

The main goal of this paper is to determine the average 
value of TZj for a large sample of disc galaxies as a function 
of their halo mass, and to use semi-analytical models to 
explore implications for how disc galaxies form. Since it is 
impossible to measure the specific angular momentum (or 
spin parameter) of a real (as opposed to simulated) dark 
matter halo, one can only obtain an estimate for the average 
value of TZj by comparing the distribution of galaxy spin 
parameters, A gal , obtained from a sample of galaxies to the 
distribution of halo spin parameters, A halo , obtained from 
numerical simulations of structure formation in a ACDM 
cosmology. This is the approach we will follow in this paper. 

Measuring A gal for an individual galaxy necessarily in- 
volves the following steps: 

(i) Determine the specific angular momentum of the 
baryons. 

(ii) Determine the mass (and virial radius) of the dark 
matter halo in which the galaxy resides. 

Step (i) is relatively straightforward. Using that, for a disc 
galaxy, 

J g ai = 2tt [ ™ E(R) V TOt (R)RdR, (4) 
Jo 

it is clear that one can determine the angular momentum of 
a disc galaxy from measurements of its surface mass density, 
E(_R), (accounting for both gas and stars), and its rotation 
curve Vrot(R) (see e.g., van den Bosch, Burkert & Swaters 
2001). Alternative, since disc galaxies in general have ex- 
ponential surface density profiles and flat rotation curves, 
one can use simple scaling relations between mass (or lu- 
minosity), size and characteristic rotation velocity in order 
to obtain an estimate of J ga i without having to measure 
E(7?) or Vrot(R) (e.g., Navarro & Steinmetz 2000; Tonini 
et al. 2006; Hernandez & Cervantes-Sodi 2006; Hernandez 
et al. 2007; Berta et al. 2008). Although less accurate than 
using Eq. Q, this method has the advantage that it can 
be applied to much larger samples of galaxies. Step (ii) is a 
little bit more involved. Typically, accurate measurements 
of halo mass require probes on scales of the virial radius. 
Examples of such 'probes' are gravitational lensing, X-ray 
data and satellite kinematics. However, unless the galaxy in 
question is the central galaxy of a cluster-sized halo, such 
measurements are rarely available for individual systems, 
and one has to rely on uncertain estimates of the halo mass 
from the galaxy's rotation curve. Alternatively, one can use 
various statistical techniques to assign to each galaxy an av- 
erage halo mass. Recent years have seen dramatic progress 
in using a variety of techniques to constrain the average rela- 
tion between galaxy properties (typically luminosity of stel- 
lar mass) and halo mass, in particular from galaxy clustering 
(e.g., Yang, Mo & van den Bosch 2003; Tinker et al. 2005), 



galaxy-galaxy lensing (e.g, Mandelbaum et al. 2006; Cac- 
ciato et al. 2009), abundance matching (e.g., Conroy, Wech- 
sler & Kravtsov 2006), and satellite kinematics (e.g., More 
et al. 2009). 

In this paper we compute the average specific angu- 
lar momentum and spin parameters of disc galaxies as a 
function of their halo mass. We improve on previous stud- 
ies in a number of ways. We constrain halo masses us- 
ing recent results for the halo mass - stellar mass relation 
for disc galaxies obtained from weak gravitational lensing 
and satellite kinematics (Dutton et al. 2010b), rather than 
abundance matching results as in Tonini et al. (2006), or 
the highly questionable assumption of a constant galaxy 
mass fraction made by Hernandez & Cervantes-Sodi (2006), 
Hernandez et al. (2007), and Berta et al. (2008). We also 
improve upon these studies by taking account of realis- 
tic uncertainties on halo masses. We consider both stel- 
lar and gaseous discs, rather than just stellar discs, which 
is important for low mass galaxies which tend to be gas 
rich (Tonini et al. 2006). We calculate disc angular mo- 
mentum by integrating a model that is constructed to re- 
produce the size-mass and velocity-mass (i.e., Tully-Fisher, 
Tully & Fisher 1977) relations, rather than simple scaling 
arguments for exponential discs in isothermal haloes (i.e., 
Jg a i/Mg a i = 2_RdKot), as used by e.g., Navarro & Steinmetz 
(2000) and Hernandez et al. (2007). We also discuss how 
uncertainties related to stellar mass-to-light ratios and adi- 
abatic contraction impact the inferred values of TZj , and use 
simple models of disc galaxy formation to discuss the im- 
plications of our findings for cooling, outflows and angular 
momentum loss. 

We adopt a flat ACDM cosmology with matter den- 
sity parameter fi m = 0.27, baryon density Qb = 0.044, and 
Hubble parameter h = Ho/lOOkms" 1 Mpc -1 = 0.7. Dark 
matter halo masses are defined so that the mean density of 
the halo (assumed to be spherical) is 200 times the critical 
density. 



2 DEFINITIONS AND MASS MODELS 

This section gives a brief overview of the main parameters 
we discuss in this paper, and the constrained mass models 
we use to constrain them. 



2.1 Mass and Angular Momentum Parameters 

There are four main parameters that are relevant for what 
follows: galaxy mass, M ga i, halo virial mass, M v j r , galaxy an- 
gular momentum, J ga i, and the angular momentum of the 
dark matter halo, J v ir. The galaxy mass is the sum of the 
mass in stars and cold gas. The virial mass is the total mass 
(stars, cold gas, hot gas, dark matter) within the virial ra- 
dius. The galaxy angular momentum is the sum of the an- 
gular momenta of the stars and of the cold gas. The virial 
angular momentum is the total angular momentum within 
the virial radius. 

Our main goal is to constrain the angular momentum 
ratio, TZj, defined as the ratio between the galaxy specific 
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angular momentum and the total specific angular momen 
turrfl 

(Jgal/Mg a l) _ A gal 



n 



(Jvir/M vir ) A^ 



(3) 



where the spin parameters A^aio and A ga i are defined by 
Eqs. {TJ and ((2J, respectively. Using that for disc galaxies 
Jgai oc -RdMgajKot, where Rd is the disc scale length and 
V TO t is a characteristic rotation velocity of the disc, one can 
use Eq. ((5]) to show that 



Rd oc TZj A! m io (t/~"~) ^ vir 



(6) 



which is the standard approach used in analytical and semi- 
analytical models to predict the sizes of disc galaxies (e.g., 
Kauffmann 1996; Mo, Mao & White 1998; Somerville & Pri- 
mack 1999; van den Bosch 2000; Cole et al. 2000). In many 
cases these studies simplify Eq. ([5]) by assuming that TZj = 1 
and that Kot/Kdr is a constant. In this paper we examine 
whether there is empirical support for the former assump- 
tion and what this implies for the formation of disc galaxies 
in general. As we will demonstrate, it is instructive to com- 
pare TZj to the galaxy formation efficiency parameter 



(M g al/M vir ) 



(7) 



(Ob/Om) 

which describes what fraction of the cosmologically available 
baryons have ended up in a galaxy. 

2.2 Mass Models 

We compute the average value for TZj for disc galaxies using 
the mass models of Dutton et al. (2011b), which have been 
constructed to reproduce the observed structural and dy- 
namical scaling relations of late-type (i.e., star-forming and 
disc-dominated) galaxies. They consist of four components: 
a stellar bulge, a stellar disc, a gas disc, and a dark matter 
halo. Both the bulge and the dark matter halo are assumed 
to be spherical and to follow a density distribution that is 
given by a Hernquist profile (Hernquist 1990) and a NFW 
profile (Navarro, Frenk, & White 1997), respectively. The 
density profile of the halo, however, is modified to account 
for the gravitational response to galaxy formation (most of- 
ten modelled as adiabatic contraction, see Blumenthal et al. 
1986). The disc is assumed to have an exponential surface 
density and to be infinitesimally thin. 

Each halo-galaxy system is described by nine parame- 
ters: three for the halo (virial mass, M v i r , concentration, c, 
and a parameter that describes the impact of galaxy forma- 
tion on the halo density profile); and six for the baryons 
(stellar mass, M sta r, gas mass, M ga s, bulge fraction, /b, 
bulge size, R^, stellar disc size, -Rd, and gas disc size, 
_R g ). The 6 parameters for the baryons are constrained by 
the observed scaling relations: Rh vs. M sta r; Rd vs. M s t ar ; 
R g /Rd vs. M s tar; /b vs. M star , and M sas /M ata r vs. M sta r (see 
Dutton et al. 2011b for details). For the relation between 
Mvir and M s tar we use the empirical relation from Dutton 
et al. (2010b), which is inferred from a large number of in- 
dependent studies and techniques, including galaxy-galaxy 

2 Note that TZj is not an efficiency because it can be larger than 
unity. 



lensing, satellite kinematics, galaxy clustering, and abun- 
dance matching. We will use the scatter among these differ- 
ent studies/techniques, as an estimate for the uncertainty in 
the Mvir vs. M s tar relation, and propagate these uncertain- 
ties to those on egf and TZj. Finally, the relation between 
halo mass, M v i r , and halo concentration, c, is taken from 
Maccio et al. (2008), and is calibrated using high-resolution 
numerical simulations. 

The usage of all these empirical and theoretical scal- 
ing relations implies that the nine model parameters can be 
reduced to two primary unknowns: the stellar initial mass 
function (IMF), and a prescription for how galaxy growth 
influences the structural properties of the dark matter halo 
(hereafter "the halo contraction model"). It is common prac- 
tice to assume that disc formation is an adiabatic process 
and that the halo responds by contracting using the adi- 
abatic contraction (hereafter AC) prescription of Blumen- 
thal et al. (1986; hereafter BFFP). However, numerical sim- 
ulations have shown that this prescription may not be suf- 
ficiently accurate, and several alternatives have been pro- 
posed (e.g., Gnedin et al. 2004; Abadi et al. 2010). We follow 
Dutton et al. (2007) and use a generalized halo contraction 
model, according to which a shell of dark matter initially 
(before disc formation) at radius ri ends up, after disc for- 
mation, at a radius 



Tf = Thpfp U . 



(8) 



Here Fbffp is the amount of contraction as predicted by the 
AC-model of BFFP, and v is a free parameter. Note that 
v = 1 corresponds to the standard adiabatic contraction 
model of BFFP, while the Gnedin et al. (2004) and Abadi 
et al. (2010) models can be well approximated by v — 0.8 
and v = 0.4, respectively. A model without contraction has 
v = 0, while v < indicates that the dark matter halo 
responds to disc formation by expanding. 

As shown in Dutton et al. (2011b), reproducing the ob- 
served Tully-Fisher relation implies different values for v for 
a different choice of the IMF. In this paper we will consider 
two different models that fit the Tully-Fisher relation (and 
all other scaling relations mentioned above) equally well: 
model A which adopts a Chabrier (2003) IMF and v = -0.5 
(i.e., halo expansion), and model B, in which we adopt the 
AC model of Gnedin et al. (2004), i.e., we set v = 0.8, and 
we assume an IMF that yields stellar mass-to-light ratios 
that are 0.3 dex lower than for the Chabrier (2003) IMF. 
As shown in Dutton et al. (2011b) both these models are in 
excellent agreement with the data. As we will show below, 
model B results in galaxy spin parameters that are ~ 30 
percent (~ 0.12 dex) higher than those of model A. This 
uncertainty is comparable to that arising from the uncer- 
tainty in the M v i r vs. M sta r relation. Throughout we will 
mainly present the results for model A, while occasionally 
showing, for comparison, the results for model B. 

For a given choice of IMF and v, the model completely 
specifies the mass distributions of disc-plus-bulge galaxies as 
a function of halo mass. From these mass models we com- 
pute the angular momentum of the baryons by evaluating 
Eq. [4] where the surface density is the sum of the surface 
densities of the gas and stars: £(-R) = E gas (-R) + £ s tar(-R). 
We assume that the bulge has no net angular momentum. 
If we assumed that bulges follow the same rotation curve as 
the disc, then this would increase the specific angular mo- 
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Figure 1. Left: Galaxy specific angular momentum vs. virial mass for disc galaxies. The shaded regions indicate the one and two 
<T uncertainties on virial mass. The dashed line shows the relation between virial specific angular momentum and virial mass from 
cosmological N-body simulations, corresponding to a spin parameter A' = 0.031. The pentagon shows estimates for the Milky Way (MW) 
with the shaded region showing one and two a uncertainties. Right: Specific angular momentum ratio (TZj) vs. virial mass. The solid 
line and shaded regions indicate the mean and the one and two a uncertainties for our fiducial model A, which assumes a Chabrier 
(2003) IMF and a halo contraction model characterized by v = —0.5. The dashed line shows the mean for model B, in which the stellar 
mass-to-light ratios are 0.3 dex smaller and v = 0.8. Note that, for both models, TZj is less than unity and approximately independent 
of virial mass. 



mentum of the galaxy by < 20%. The effect is largest for the 
most massive galaxies, but since massive bulges, in general, 
are not rotationally supported, we believe our assumption 
that bulges have no angular momentum to be of no major 
concern for what follows. 

A shown by Tonini et al. (2006) a significant fraction of 
the angular momentum of a galaxy is in the cold gas. In our 
models the cold gas in high mass galaxies (M ga i ~ 10 11 M Q ) 
contains ~ 30% of the angular momentum. This fraction 
rises to over 50% for galaxy masses below M ga i ~ 1O 1O M0. 
Thus if cold gas fractions increase at higher redshifts, as is 
commonly thought (but see Dutton et al. 2010a), measur- 
ing angular momentum of high redshift galaxies will require 
measurements of the distribution of cold gas, and not just 
the stars. 



3 RESULTS 

Using our constrained mass models we now discuss a number 
of correlations between mass and angular momentum of disc 
galaxies and dark matter haloes, and compare with previous 
results. These correlations are presented in Figs. IU3I and 
Tabled] 

3.1 Galaxy Mass and Specific Angular 
Momentum 

The left-hand panel of Fig. ffj shows the correlation be- 
tween galaxy specific angular momentum and virial mass 
for model A. The shaded regions show the one and two a 



uncertainties in virial masses at fixed galaxy mas^f], which is 
the dominant source of uncertainty. The dashed line shows 
the relation between virial specific angular momentum and 
virial mass from cosmological N-body simulations (Maccio 
et al. 2008), which corresponds to a spin parameter of 
A' = 0.031. The galaxy specific angular momentum scales 
with virial mass as J ga i /M ga i oc M°;J° , which is very similar 
to the virial scaling of J v \t / M V - 1T oc M^-'f 7 . Thus the ratio TZj 
between the galaxy and virial specific angular momentum is 
roughly constant (right-hand panel of Fig. [TJ. However, the 
galaxy specific angular momentum is significantly lower than 
the virial specific angular momentum with TZj = 0.611q 11 
(la) for a virial mass of M v i r = 10 Mq. The dashed line in 
the right-hand panel of Fig. [T] corresponds to model B, and 
is shown for comparison. Note that the combination of halo 
contraction and lower stellar mass-to- light ratios result in TZj 
values that are ~ 0.1 higher than for our fiducial model A. 
This difference is comparable to the la uncertainty on TZj 
due to the uncertainties in the relation between stellar mass 
and halo mass. 

The left-hand panel of Fig. [2] shows the correlation be- 
tween galaxy mass and virial mass for the disc galaxies con- 
sidered here. The dashed line shows the relation for the cos- 
mically available baryonic mass, i.e., M ga i = (f2b/fl m )-?Wvu-, 
corresponding to egf = 1. The right-hand panel of Fig. [2] 
shows the galaxy formation efficiency, egf, which is the 
ratio between the solid and dashed lines in the left-hand 

3 Note that these are the uncertainties on the mean; they do not 
reflect the scatter. 
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M*/[M.] M vir /[MJ 

Figure 2. Left: Galaxy mass vs. virial mass for disc galaxies. The shaded regions reflect the one and two <r uncertainties on virial masses. 
The dashed line corresponds to the universal baryon fraction (f2b/fi m = 0.163). The pentagon shows estimates for the Milky Way (MW), 
with the shaded region showing one and two a uncertainties. Right: Galaxy formation efficiency (cgf) vs - virial mass. As in Fig. [T] wc 
show the results for both models A and B. For clarity, we only show the impact of virial mass uncertainties for our fiducial model A. 
Note that unlike IZj , which reveals no significant dependence on halo mass (see Fig. [TJ , the galaxy formation efficiency varies by a factor 
of ~ 3 over the mass range probed. 



panel. Unlike the angular momentum ratio, the galaxy for- 
mation efficiency is a strong function of virial mass, varying 
by a factor of ~ 3 from egf — 0.1 at the low mass end 
(M vir ~2x lO n M ) to e G F ^ 0.3 at the high mass end 
(M vir ~5x 10 12 M Q ). Once again, the dashed line in the 
right-hand panel corresponds to model B, which yields lower 
galaxy formation efficiencies due to the lighter IMF. 

As we demonstrate and discuss below, (semi)-analytical 
models and hydrodynamical simulations of galaxy formation 
do generally not predict values of IZj that are independent 
of halo mass, and if they do, the corresponding values of egf 
are inconsistent with the empirical results inferred here. 

3.2 Galaxy Spin Parameter 

We now discuss the galaxy spin parameter, as this is a quan- 
tity that has been measured and interpreted by previous 
studies. Fig. [3] shows the average galaxy spin parameter as 
a function of halo virial mass. This is the same as the right- 
hand panel of Fig. [1] but with a zero point shift. As discussed 
above, there is a slight mass dependence to A gal , but given 
the uncertainties in halo masses, A gal is consistent with be- 
ing independent of halo mass. The normalization is ~ 60% 
lower than predicted for ACDM haloes (dashed line). There 
are two simple interpretations of this result. One is that the 
baryons that form disc galaxies have lost about 40% of their 
angular momentum during galaxy formation, which could 
occur via dynamical friction. Another is that the baryons 
that form disc galaxies have 40% lower specific angular mo- 
mentum than the cosmically available baryons, which could 
occur due to feedback and/or inefficient cooling. We discuss 
these interpretations in more detail in § [4j 



An interpretation advocated by Tonini et al. (2006) is 
that if one restricts attention to haloes that have not expe- 
rienced a major merger since z ~ 3, the average spin param- 
eters A' turns out to be around 0.023 (D'Onghia & Burkert 
2004). This is motivated by the idea that bulge- less disc 
galaxies require a quiescent merger history. However, this 
simple expectation has been cast into doubt by numerical 
simulations which have shown that (gas rich) discs can sur- 
vive major mergers (e.g., Springel & Hernquist 2005; Hop- 
kins et al. 2009). While we cannot rule out a merger history 
bias for massive disc galaxies, we argue that such a bias 
cannot occur for low mass galaxies (M ga i < 3 x 1O 1O M0, 
M v i T ^ 1O 12 M0) because typical low mass galaxies are 
bulge-less, and thus must form in haloes with typical merger 
histories (and hence typical halo spins). 



3.3 Comparison with the Milky Way 

As a reference, the black pentagons and light blue shaded 
regions in Figs. [T] & [5] show estimates for the specific an- 
gular momentum and mass of the Milky Way (MW). Here 
we have estimated the specific angular momentum of the 
MW assuming J ga i/Af ga i = 2i?dKot, which is appropriate 
for an exponential disc with a flat rotation curve. We adopt 
Rd = 2.8 ± 0.23kpc(lo") based on the dynamical models 
of Widrow et al. (2008) and V ro t = 219 ± 20 kms^lcr) 
from the observations of Reid et al. (1999). For the baryonic 
mass of the MW we use the bulge and disc masses, again 
from the dynamical models of Widrow et al. (2008) to obtain 
log 10 (M ga i/M Q ) = 10.71T0.075 (la). For the virial mass of 
the MW we average the results of Smith et al. (2007) and 
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Table 1. Relations between mass and angular momentum of galaxies and dark matter haloes for model A, as presented in Figs, Hl3l 
Uncertainties on halo masses are 1 and 2a, and are propagated to other parameters as appropriate. 
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Figure 3. Average spin parameter vs. halo virial mass, for disc 
galaxies (solid line and grey shaded region), and ACDM dark mat- 
ter haloes (dashed line, red shaded region). For the galaxy spin 
parameter the shaded regions show the one and two a uncertain- 
ties introduced by the systematic uncertainty in halo masses. The 
pentagon shows an estimate for the Milky Way (MW) with the 
shaded region showing one and two a uncertainties. 



Xue et al. (2008) yielding log 10 (M v ir/M Q ) = 11.97 ± 0.22 
(2a). 

The resulting galaxy spin parameter (see Eq.2) and 
galaxy formation efficiency (see Eq.7) of the MW are Ag al = 
0.03ll^oii (2cr), and e GF = 0.033±g:?| {2a), respectively. 
Both Ag al and egf are higher than for typical late- type galax- 
ies. This could signify that the MW is atypical, or it could 
simply be a reflection of the measurement errors (especially 
in halo masses) both for the MW and for samples of late- 
type galaxies. 



3.4 Comparison with Previous Studies 

Several studies in the past used similar methods to deter- 
mine the spin parameters of disc galaxies. These studies all 
find higher spin parameters than we find here, but as we 
show the differences can be ascribed to differences in the 
estimation of halo masses, which is the dominant source of 
systematic uncertainty in measuring the spin parameter. 

Van den Bosch, Burkert & Swaters (2001; hereafter 
BBS01) used Eq. (Q and halo masses inferred from rota- 
tion curve modelling to determine the galaxy spin parame- 
ters for a sample of 14 disc-dominated dwarf galaxies. They 
found that (A gal ) = 0.048 (see Fig. gfl, albeit with a rel- 
atively large uncertainty due to the small sample size and 
the uncertainties in the stellar mass-to-light ratios and halo 
masses. At first sight, this is much larger than the average 
galaxy spin parameters inferred here; (A gal ) = 0.019. How- 
ever, there are several factors that can account for this large 
difference. BBS01 used dark halo masses and sizes derived 
from fits to rotation curves assuming adiabatic contraction. 
Since the rotation curves only extend to ~ 10 — 20% of the 
halo virial radius, this involves a significant extrapolation, 
and thus the potential for systematic errors. As discussed 
in detail in Dutton et al. (2007), the assumption that disc 
formation results in adiabatic contraction of its dark matter 
halo may not be realistic. The galaxies in BBS01 have a me- 
dian Vi as t/Vvir = 1.55, where VJ ast is the rotation velocity at 
the last measured point (i.e., at the largest radius) . Since 
dark matter haloes in the ACDM concordance cosmology 
have Knax/Vvir J; 1-2, where V m ax is the maximum circular 
velocity, the derived values from BBS01 are strongly influ- 
enced by the assumption of adiabatic contraction. Dutton 
et al. (2010b) have shown that the halo virial velocity is, on 
average, approximately equal to the optical rotation veloc- 
ity. Using this assumption to derive halo masses results in 
log lo (M v i r /M0) = 11.2 ± 0.3 and a median spin parameter 
A gal = 0.023 (right-hand panel of Fig. [4| . If we consider a 
10% uncertainty on V I0 t/V v i T (Dutton et al. 2010b) then this 



4 Here we have converted the results from BBS01 to the spin 
parameter definition of Bullock et al. (2001). 
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Figure 4. Histograms (hatched) of the distribution of galaxy spin 
parameters (Ag al ) f° r the 14 disc dominated dwarf galaxies from 
van den Bosch, Burkert, & Swaters (2001; hereafter BBS01). In 
the left-hand panel A gal has been obtained using the halo virial 
velocities from BBS01. In the right-hand panel, the halo virial ve- 
locities have been calculated assuming V vlr = Vi as t , where \\ aat is 
the rotation velocity at the last measured point (i.e., at the largest 
radius). The latter results in spin parameters in good agreement 
with our results (cf. Fig.[3jl- 
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results in a range of 0.019 — 0.028 for A ga i, which is consistent 
with our measured values (see Fig. |3J). 

Tonini et al. (2006) used average scaling relations be- 
tween stellar mass, gas mass, halo mass and disc size to 
compute an average galaxy spin parameter as function of 
halo mass, under the assumption that dark matter haloes 
follow a Burkert (1995) profile. They found that (A gal ) ~ 
0.025 — 0.030 with no significant dependence on halo mass. 
A comparison between the halo mass - stellar mass relation 
used by Tonini et al. (2006) and by ourselves indicates that 
their halo masses are a factor of ~ 2 lower than ours. This 
difference in halo masses fully accounts for their higher spin 
parameters. 

Finally, Hernandez et al. (2007) determined the galaxy 
spin parameters for a large sample of galaxies taken from 
the Sloan Digital Sky Survey (SDSS). Rather than Tonini 
et al. (2006), they obtained estimates for A ga i for individ- 
ual galaxies using the method proposed by Hernandez & 
Cervantes-Sodi (2006), which is based on the following as- 
sumptions: (i) dark matter haloes are singular isothermal 
spheres, (ii) the self-gravity of the disc can be ignored, and 
(iii) all galaxies have the same ratio of dark matter mass to 
disc mass. In addition, they inferred the rotation velocity of 
each individual galaxy in a statistical sense, using an empir- 
ical Tully-Fisher relation. They found that for disc galaxies 
A gal follows a log-normal distribution with (A gal ) ~ 0.059, 
which is a factor of ~ 3 higher than our result. The method 
used by Hernandez et al. (2007) makes several highly ques- 
tionable assumptions, and thus we do not consider their 
result reliable. For example, Hernandez et al. (2007) as- 
sume a constant galaxy mass to dark matter mass ratio of 
F = 0.04, which corresponds to a galaxy formation efficiency 
of egf — 0.25. Comparison with Fig. [2]shows that while this 
is a good assumption for a 2 x 1Q 12 Mq halo, it significantly 
over-estimates the galaxy formation efficiency, and hence 
spin parameters, of lower mass galaxies. Since Hernandez 
et al. (2007) use a volume limited sample, their measure- 
ment of the spin parameter is biased towards lower mass 
galaxies, where their assumption of F = 0.04 is high by a 
factor of ~ 2. In the method used by Hernandez et al. (2007) 



Figure 5. The cumulative, normalized fractions of mass (blue 
line) and angular momentum (dashed red line) for ACDM haloes 
(Navarro et al. 1997; Bullock et al. 2001). The shaded region 
shows the uncertainty on J(r)/J v ; r assuming s = 1.3 ± 0.3 in 
Eq. This clearly demonstrates that mass is more centrally 
concentrated than angular momentum. For example, the inner 
50% of the mass contains just 20_t^% of the angular momentum 
(black dotted lines). Thus if the efficiency of galaxy formation is 
regulated by the efficiency of cooling, we would expect a strong 
correlation between galaxy spin parameter and galaxy formation 
efficiency, which is not observed. 



A gal oc F, and thus their assumption of F = 0.04 accounts 
for at least half the discrepancy between our result and 
theirs. 



4 COMPARISON WITH GALAXY 
FORMATION MODELS 

4.1 Inside-Out Cooling Models 

If the galaxy formation efficiency is determined by the ef- 
ficiency of cooling, then we expect that IZj will be higher 
for higher galaxy formation efficiencies. This is because in 
ACDM haloes mass is more centrally concentrated than an- 
gular momentum. This is illustrated in Fig. [5] The blue solid 
line shows the cumulative mass fraction, M(< r)/M v i t as a 
function of radius for a NFW halo with concentration c = 10. 
The red dashed line shows the cumulative angular momen- 
tum fraction, J(< r)/J v i r , as a function of radius assuming 



J(< r)/J vir 



[M(< r)/M vir ] fl 



(9) 



M(< r)/M vl 

where s = 1.3±0.3 (Bullock et al. 2001). Under the assump- 
tions that (i) the baryons initially acquire the same specific 
angular momentum as the dark matter, and (ii) the baryons 
accrete onto the central galaxy inside-out, then Eq.[9]implies 

T^i ~ e GF- 

We now extend this simple inside-out cooling model to 
cosmologically evolving dark matter haloes. The red lines in 
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Figure 6. Correlations between halo mass, M v i r , galaxy forma- 
tion efficiency, eQFi an d specific angular momentum ratio, IZj, for 
models and observations at redshift z = 0. The empirical relations 
for disc galaxies are shown as solid black lines, with the shading 
reflecting the one and two a uncertainties on halo masses. Note 
that the galaxy formation efficiency increases with increasing halo 
mass, while the specific angular momentum ratio is independent 
of halo mass. The corresponding relations for a simple cooling 
model are shown with red lines, where the different line types 
correspond to different metallicities, as indicated. In this cooling 
model IZj ~ £qf (independent of the metallicity of the hot gas), 
inconsistent with the empirical findings. 

Fig.[6]show the relations between IZj , egf and M v \t for a sce- 
nario in which the galaxy formation efficiency is determined 
by the efficiency of cooling. This model is a simplified version 
of the galaxy formation model from Dutton & van den Bosch 
(2009) . Briefly this model consists of NFW haloes that grow 
by smooth accretion of baryons and dark matter. The angu- 
lar momentum of each shell of accreted gas is determined by 
assuming that the angular momentum distribution (AMD) 
of a given halo is independent of time. The AMD is spec- 
ified by two parameters: the normalization (A) and shape 
(a). When the baryons enter the halo they are shock heated 
to the virial temperature. The cooling time depends on the 
temperature, density and metallicity of the hot gas. We use 
the metallicity dependent collisional ionization equilibrium 
cooling functions of Sutherland & Dopita (1993). Results for 
four different metallicities from primordial to one third so- 
lar are shown. For low mass haloes (M v i r ~ 10 11 Mq) cool- 
ing is very efficient, so that ~ 90% of the baryonic mass 
is accreted (left-hand panelsjf). For more massive haloes, 
the galaxy formation efficiency, egf, depends strongly on 
the metallicity of the gas, with higher metallicities yielding 
larger values of egf- Clearly, the trend that egf decreases 
with increasing halo mass is completely opposite to the em- 



5 It is not 100% because it takes a free fall time for the most 
recently accreted baryons to reach the galaxy. 



Figure 7. Same as Fig. [6] but this time we compare the empirical 
relations (solid black lines with shaded regions) with models that 
include feedback and/or angular momentum loss (red lines). All 
models shown adopt a hot gas metallicity of Z-^ot = O.IZq, but 
they differ in feedback efficiency, epBi an d the fraction of angular 
momentum lost by the baryons during galaxy formation, /i oa t. 
The dotted line (epe = /i os t = 0) is identical to the simple cooling 
model from Fig. [6] and is shown for comparison. Note that none 
of the models is able to simultaneously match the low galaxy 
formation efficiencies and the constant angular momentum ratio. 
See text for a detailed discussion. 

pirical result. As for the parameter TZj, the simple cooling 
models predict that TZj should decrease with increasing halo 
mass, which is a consequence of cooling being an inside-out 
process combined with the fact that mass is more centrally 
concentrated than angular momentum (cf. Fig. [5| . Hence, in 
massive haloes where egf -C 1 due to inefficient cooling, the 
galaxy ends up having less specific angular momentum than 
its dark matter halo. This strong mass dependence of IZj 
is again inconsistent with our empirical results. Finally, the 
cooling models predicts an almost linear relation between IZj 
and £gf, which is independent of the gas metallicity (lower 
right-hand panel of Fig. [6]). Clearly, this linear relation is 
very different from its empirical equivalent, indicating that 
additional mechanisms for accreting and/or expelling gas 
are required. 

4.2 Outflows and Angular Momentum Loss 

We now consider two modifications to the simple cooling 
model: outflows and angular momentum loss. We assume 
that outflows, which move at the local escape velocity, are 
driven by energy from supernovae (SNe). The fraction of 
the SNe energy that goes into an outflow is specified by 
the free parameter epe (see van den Bosch 2002; Dutton 
& van den Bosch 2009). Other than inefficient cooling in 
massive haloes, and UV heating in low mass haloes, galac- 
tic outflows are the primary mechanism for explaining the 
observed low galaxy formation efficiencies (e.g., Dekel & 
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Silk 1986; Cole et al. 1994). For simplicity, we assume that 
outflow gas leaves the disc and halo, and does not return. 
However, recent cosmological hydrodynamical simulations 
have shown that ejected gas may return to the galaxy (Op- 
penheimer et al. 2010). Some of the ejected gas may mix 
with the halo gas and thus re-accrete with higher specific 
angular momentum (Brook et al. 2012). Since this process 
is a re-distribution (i.e., ejected gas gains angular momen- 
tum, halo gas loses angular momentum), it does not change 
the specific angular momentum of the baryons. If the halo 
gas accretes onto the central galaxy (as is expected for low 
mass, M v i r < 1O 12 M0, haloes) then there is no net change 
in the angular momentum of the galaxy. If the halo gas does 
not accrete (which is expected for high mass haloes) then 
the galaxy can gain angular momentum via this galactic 
fountain effect. Thus if we were to implement re-accretion 
into our model we would need higher feedback efficiencies to 
achieve the same galaxy formation efficiencies, but to first 
order there would be no change to the specific angular mo- 
mentum. 

We assume that during the accretion process, a frac- 
tion /i ost of the angular momentum of the gas is lost to the 
dark matter halo. The mechanism responsible for this trans- 
fer is dynamical friction. In our models we do not take into 
account the effects of angular momentum transfer and out- 
flows on the structure of the dark matter halo. But we note 
that both of these processes will result in halo expansion, or 
at least a reduction of the amount of adiabatic contraction, 
i.e., a lower effective value of v (e.g., Navarro et al. 1996; El- 
Zant et al. 2001; Mo & Mao 2004; Cole, Dehnen & Wilkinson 
2011) 

In Fig. we show several models with different values 
of epB and /i os t- In all models we assume that halo gas is 
enriched to one tenth solar. The dotted line shows a model 
with efb = 0.0 and /i os t = 0.0: this corresponds to the sim- 
ple cooling model from Fig. [5] and is shown for comparison. 
The galaxy formation efficiencies are too high for all halo 
masses. In order to produce realistic galaxy formation effi- 
ciencies outflows are needed. A model with efb = 0.5 and 
/lost = 0.0 (short dashed line) has egf ~ 20%. In our model 
outflows increase IZj by preferentially removing gas with 
low specific angular momentum. Since outflows are driven 
by star formation, regions that have higher star formation 
rates have higher mass outflow rates. Star formation is em- 
pirically more efficient at smaller galactic radii (Kennicutt 
1998) as well as at higher redshifts (e.g., Noeske et al. 2007; 
Daddi et al. 2007; Elbaz et al. 2007). Smaller galactic radii 
correspond to lower specific angular momentum and discs 
are smaller at higher redshifts (e.g., Dutton et al. 2011a); 
both of these effects contribute to outflows increasing the 
specific angular momentum of the material that remains in 
the galaxy. In our energy driven outflow model the mass 
loading factor (= outflow rate / star formation rate) is higher 
in lower mass galaxies, which results in IZj increasing more 
in lower mass haloes (by up to a factor two). This increases 
the disagreement with the empirically inferred relation be- 
tween IZj and halo mass. We can lower IZj by introducing 
angular momentum loss (i.e., setting /i ost > 0). The long 
dashed-dotted lines in Fig. [7] correspond to a model with 
epB = 0.5 and /i ost = 0.5. This angular momentum loss re- 
duces IZj by a factor ~ 2 at all mass scales, bringing it in rea- 
sonable agreement with the observed value for haloes with 
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Figure 8. Same as Fig. [7] but this time we have tuned the feed- 
back efficiency, epe, and angular momentum loss, /i ost , of the 
model (red, solid dots) in order to accurately reproduce the em- 
pirical relations (solid black lines with shaded regions). The small, 
upper right-hand panel shows the corresponding halo mass depen- 
dence of epB (solid circles) and /i ost (open circles): reproducing 
the empirical relations requires both of these parameters to de- 
crease strongly with increasing halo mass. See text for a detailed 
discussion. 
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Figure 9. Same as Fig. [8] but here we have used model B (stellar 
mass-to-light ratios that are 0.3 dex lower and halo contraction 
characterized by v = 0.8) for the empirical relations. Reproduc- 
ing the inferred galaxy formation efficiencies now requires even 
stronger feedback (i.e., larger £fb)i while the fraction of angular 
momentum that is lost is similar to that for our fiducial model 
shown in Fig. [8] 
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M v i r ~ 1O 12 M0. However, for lower (higher) mass haloes 
the model overpredicts (underpredicts) the value of TZj com- 
pared to the empirically inferred values. Increasing /i oat also 
modified the galaxy formation efficiency because it results in 
denser (smaller) discs. In low mass haloes, the galaxies are 
gas rich. Hence, an increase in disc density results in higher 
star formation rates and consequently more energy to drive 
outflows. This more than compensates for the higher escape 
velocities, resulting in lower values of 6gf- In massive haloes, 
on the other hand, the galaxies are gas poor, so that an in- 
crease of the density of the disc only results in an increase 
of the escape velocity. This reduces the amount of ejected 
gas, and therefore increases egf- 

It is clear that none of the models discussed thus far 
can simultaneously reproduce the empirically inferred trends 
of £gf and TZj with halo mass. As an illustration of what 
seems to be required, we now consider a model in which both 
£fb and /i os t depend on halo mass. The red filled symbols 
in Fig. [8] show a model in which efb and /i ost have been 
tuned so that the model matches the empirical cgf and TZj . 
The corresponding values of efb and /i os t, as a function of 
virial mass, are shown in the upper right-hand panel. Lower 
mass haloes require higher efb and /i os t. For example, this 
model has (e F B,/iost) = (0.7,0.7) for Af vir ~ 10 113 M Q , and 
(e F B,/lost) = (0.25,0.0) for M vir ~ 10 12A M e . 

The results shown in Fig. [8] correspond to the empirical 
£gf and TZj inferred under the assumptions of model A (i.e., 
Chabrier IMF and halo contraction characterized by v — 
—0.5). Fig.[9]shows the analogous results but for model B (an 
IMF that results in stellar mass-to-light ratios that are 0.3 
dex lower than for model A and with v = 0.8). Reproducing 
the galaxy formation efficiencies now requires even stronger 
feedback, but with a similar mass dependence, whereas the 
fraction of angular momentum that is lost is similar to that 
of our fiducial model. 

4.3 Disc Stability 

Thus far we have assumed that disc galaxies form in typical 
dark matter haloes. This must be true at the low mass end, 
simply because disc galaxies dominate the galaxy population 
at low masses. However, at high masses disc galaxies are a 
minority, so that it is possible that they form in a biased 
subset of dark matter haloes. To reconcile our models with 
the empirical fact that TZj seems to be independent of halo 
mass requires a bias such that massive disc galaxies form 
preferentially in haloes with larger spin parameters. 

Such a bias is naturally achieved by invoking disc stabil- 
ity. Disc galaxies can only survive as such to the present day 
if their discs are sufficiently stable. Disc galaxies that are 
unstable redistribute their angular momentum via secular 
evolution, typically resulting in spheroid dominated galax- 
ies (e.g., Combes et al. 1990; Norman, Sellwood & Hasan 
1996; van den Bosch 1998; Mao & Mo 1998). Since more 
self-gravitating (denser) discs are more unstable (e.g., Ef- 
stathiou, Lake & Negroponte 1982), and since the surface 
density of a disc is proportional to A^ lo (e.g., Mo, Mao 
& White 1998), one expects a bias against disc galaxies in 
haloes with the lowest spin parameter. Since the disc's sur- 
face density is also proportional to 6gf , and empirically egf 
increases with increasing halo mass, one expects discs in 
less massive haloes to be more stable than those in more 




10" 10' 2 10" 10" 10" 10" 

Figure 10. The impact of disc stability on the galaxy forma- 
tion efficiency (cqf; upper left-hand panel) and specific angu- 
lar momentum ratio (TZj; lower let-hand panel). The solid line 
with the shaded regions indicate the empirical relations (for our 
fiducial model A). The dashed, red line corresponds to a simple 
model with Zi lot = O.IZq, epB = 0.5 and /i ost = 0.6, which 
yields galaxy formation efficiencies in reasonable agreement with 
the data. However, it predicts that TZj decreases strongly with 
increasing halo mass, in clear conflict with the data. The dot- 
dashed, black line corresponds to the same model, but now we 
have removed those discs that do not satisfy our disc stability 
criterion. This removes massive galaxies in low-spin haloes, re- 
sulting a higher TZj for mass haloes, in better agreement with the 
data. 



massive haloes. Hence, disc stability may introduce a bias 
against forming massive disc galaxies in low spin parameter 
haloes, exactly what seems to be required. 

To test the potential impact of disc stability on the 
relations between sof, TZj, and M v j r , we add a simple sta- 
bility criterion to our models: we assume that if the disc 
contributes more than 84% of the mass within 2.2 disc scale 
lengths it will be highly unstable, and we remove these galax- 
ies from our sample of model galaxies. We note that discs 
are already expected to become unstable — and hence to 
form bulges via secular evolution — at disc mass fractions 
lower than 84%. However, we choose a high stability thresh- 
old as we do not wish to remove the galaxies that remain 
disc dominated after secular evolution. 

Fig. [10] shows the effect of disc stability on egf and 
TZj. In our fiducial model (red long-dashed lines) the frac- 
tion of gas that cools at each time step (/ coo i) depends on 
the metallicity of the hot halo gas (here Zhot = O.IZq), 
the feedback efficiency (here £fb = 0.5) and the fraction 
of angular momentum that is lost (here fy ost — 0.6). This 
model approximately reproduces the observed relation be- 
tween galaxy formation efficiency and halo mass, but it pre- 
dicts that TZj decreases strongly with increasing halo mass, 
which is inconsistent with the empirical results. Invoking the 
disc stability criteria discussed above (black dot-long-dashed 
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lines) results in a lower fraction of stable discs in higher mass 
haloes (lower right panel): when computing the average TZj 
at a given halo mass (for both models and observations), 
we assume that the haloes have a median spin parameter 
as expected from cosmological simulations (i.e., A' = 0.03f). 
However, since stable discs have higher median spin param- 
eters than unstable discs, this results in stable discs having 
higher TZj. The stable discs in our model still have an anti- 
correlation between TZj and halo mass, but the model is now 
consistent with the observations at the 2a level. 

The disc stability scenario also provides a qualitative 
explanation for two additional observational facts: (i) the 
scatter in disc sizes decreases with increasing stellar mass 
(Shen et al. 2003), such that for the most massive discs the 
scatter is half that expected from the scatter in halo spin 
parameters (Dutton et al. 2011a); and (ii) the fraction of 
galaxies that are star forming disc galaxies declines with 
increasing galaxy and halo mass (e.g., Yang et al. 2008). 
Hence, we conclude that disc stability may be an important 
ingredient for understanding the empirical scaling relations 
between egf, TZj and halo mass. 
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4.4 Comparison with Hydrodynamical 
Simulations 

The simplistic analytical models used above indicate that it 
is not easy to understand why the specific angular momen- 
tum ratio TZj seems to be independent of halo mass, whereas 
the galaxy formation efficiency increases strongly with halo 
mass. In fact, reproducing the empirical scaling relations be- 
tween £gf, TZj and halo mass seems to require a feedback 
efficiency that declines with increasing halo mass. 

An important downside of our simplistic models is that 
they do not properly account for the hydrodynamics of out- 
flows from galaxies embedded in a large scale environment 
from which the galaxy is also accreting matter. This requires 
cosmological, hydrodynamical simulations, which arguably 
are the best tool available to model the complex, hierarchi- 
cal nature of galaxy formation. Recently, Sales et al. (2009) 
have presented a suite of high resolution, hydrodynamical 
simulations of the formation of galaxies in a ACDM con- 
cordance cosmology, with a wide range of feedback recipes. 
Their resulting galaxies (at z = 2) have properties that are 
well fit by 

TZj = 9.71m ga i[l - exp(-l/9.71m ga i)]. (10) 

independent of the feedback recipe! Here m ga i = 
Mg a i/M v i r = £gf (fib/^m) is the galaxy mass fraction. As a 
specific example, for a galaxy mass fraction of m ga i = 0.033 
(i.e. 20% of the cosmic baryon fraction) the specific angular 
momentum ratio is predicted to be TZj = 0.31. This scaling 
is given by the black dashed line in Fig. 1111 Interestingly, the 
relation from Sales et al. (2009) is very similar to the results 
from our simple cooling model (red lines) . This suggests that 
feedback in their simulations largely preserves the rank or- 
der in binding energy of the baryons, such that the relation 
between egf and TZj of the resulting galaxies is predomi- 
nantly set by the relative radial profiles of mass and angular 
momentum shown in Fig. [5] It is unclear at present why the 
simulations seem to preserve the rank order of binding en- 
ergy, but it may have to do with the fact that outflows are 



Figure 11. Correlation between specific angular momentum ra- 
tio, TZj, and galaxy formation efficiency, egp. The red lines show 
the relations from our simple inside-out cooling model with dif- 
ferent gas metallicities (see Fig.[6]l. The black, dashed line shows 
the relation from Sales et al. (2009), obtained for simulated galax- 
ies at z ~ 2 in a cosmological hydrodynamical simulation. Note 
that these simulated galaxies follow exactly the relation between 
TZj and eqf expected for a simple cooling model, despite the fact 
that the simulations include supernova feedback. See text for a 
detailed discussion. 

inhibited from travelling far from their source due to pres- 
sure confinement of the surrounding gas, an effect that is 
not accounted for in our simplistic, analytical models. Irre- 
spective of what causes the relation between TZj and egf in 
the simulations, it is clear that Eq. [10] is inconsistent with 
the empirical relation. We caution, though, that the simula- 
tions were only run to redshift z = 2, whereas our empirical 
results correspond to z ~ 0. In addition, our empirical re- 
sults are for star-forming disc-dominated galaxies, whereas 
the results from the simulations are for all types of galaxies. 
It remains to be seen whether simulations of disc formation 
that are run to z = yield a relation between TZj and egf 
in better agreement with the data. 



5 SUMMARY 

We have combined measurements of halo virial masses from 
weak lensing and satellite kinematics with the Tully-Fisher 
and size- mass relations of disc galaxies (i.e., star- forming 
and disc-dominated) to infer the average spin parameter of 
disc galaxies as a function of their halo mass. Using toy mod- 
els for disc galaxy formation we have discussed implications 
of these results for cooling, outflows and angular momentum 
loss. We summarize our results as follows: 

• The average galaxy spin parameters of disc galaxies are 
consistent with being independent of halo mass for the range 
of halo masses probed here: 11.3 ^ log 10 (M v i r /-/Wo) 12.7. 
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• The primary uncertainty in measuring the spin param- 
eters of galaxies is the determination of halo masses. Real- 
istic uncertainties of 0.25 dex (2a), result in uncertainties 
in galaxy spin parameters of a factor of 1.5. A secondary 
uncertainty is the stellar IMF. Since gas discs have higher 
specific angular momentum than stellar discs, lower stellar 
mass normalizations will result in higher spin parameters. 
For example, stellar masses lower by 0.3 dex result in spin 
parameters higher by 0.12 dex. 

• The average spin parameter of galaxies, (A ga i) = 
0.019tQQQ3 (la), is smaller than that of their host dark 
matter haloes, (Ahalo) = 0.031 ± 0.001 (Maccio et al. 2007, 
2008; Bett et al. 2007). The inferred spin parameters of disc 
galaxies reveal no dependence on halo mass, so that the 
specific angular momentum ratio TZj ~ 0.6, independent of 
halo mass. The galaxy formation efficiency parameter, egf, 
however, reveals a strong mass dependence, increasing from 
~ 10% for M vir = lO n - 5 M to ~ 30% for M vir = 10 12 - 5 M© 
for our fiducial model. 

• Since mass is more centrally concentrated than angu- 
lar momentum in ACDM haloes a simple inside-out accre- 
tion model results in a strong correlation between galaxy 
spin parameter and galaxy formation efficiency, contrary to 
observations. This provides further support for the already 
popular notion than the low galaxy formation efficiencies in 
haloes of mass M V - 1T < lO n ' 7 M0 are determined primarily 
by feedback processes, and not by an inefficiency of cooling. 

• The empirically inferred relations between galaxy for- 
mation efficiency, 6gf, angular momentum ratio, TZj, and 
halo mass seem to indicate that galaxy formation involves 
galactic outflows, or other mechanisms that can regulate the 
galaxy formation efficiency (e.g., supernova feedback, AGN 
feedback, reionization, pre-heating) , as well as mechanisms 
that cause a substantial transfer of angular momentum from 
the baryons to the dark matter (e.g., dynamical friction). 
Although the need for feedback mechanisms to regulate the 
efficiency of galaxy formation has long been recognized, the 
observed scaling relations for disc galaxies seem to require 
that the often adopted supernovae feedback efficiency pa- 
rameter, £fb, decreases with halo mass. In addition, the 
angular momentum content of observed disc galaxies also 
require that angular momentum loss is more important for 
disc formation in less massive haloes. It remains to be seen 
whether realistic models for galaxy formation can achieve ef- 
fective values for efb and /i os t that reveal such a dependence 
on halo mass. 

What could be the cause for the effective values for efb 
and /i ost to decrease with increasing halo mass? One possi- 
bility is that this mass dependence is a manifestation of the 
transition from cold-mode accretion at M vlv < lO n ' 7 M0 to 
hot-mode accretion for haloes with M v \r ^ lO n ' 7 M0 (e.g., 
Birnboim & Dekel 2003; Keres et al. 2005, 2009; Brooks 
et al. 2009). In the cold mode regime, outflows are unhin- 
dered by hot-gas atmospheres, which may result in higher 
effective feedback efficiencies. If, as envisioned in our models, 
the outflows preferentially remove low angular momentum 
material, the same outflows which reduce the galaxy forma- 
tion efficiency are expected to increase the specific angular 
momentum ratio, TZj. Hence, in order to reproduce the in- 
ferred IZj ~ 0.6, the baryons that end up in the disc need to 
lose a significant fraction (~ 60%) of their angular momen- 



tum during the galaxy formation process. If the cold accre- 
tion is sufficiently clumpy, this may come about because of 
dynamical friction. 

In the hot mode accretion regime, outflows are expected 
to be less efficient, since they have to do work against their 
hot gaseous atmospheres. In this regime, one expects that 
the baryon fraction and angular momentum of the resulting 
disc galaxies are largely determined by the inside-out cool- 
ing of the hot gas, resulting in IZj decreasing strongly with 
increasing halo mass. This, however, is inconsistent with the 
empirical result that IZj appears to be independent of halo 
mass. We have argued that disc stability may play an impor- 
tant role here, effectively removing systems with low halo 
spin parameter from the sample of disc galaxies (because 
they produce unstable discs), thereby increasing the aver- 
age IZj of massive disc galaxies. Another important process 
may be that at high redshifts (z > 2), massive haloes in 
the hot-mode accretion regime can still accrete a significant 
amount of material via cold streams that penetrate the hot 
halo (e.g., Dekel et al. 2009). Depending on the 'impact pa- 
rameters' of such streams, the cold material thus deposited 
to the central galaxy may have specific angular momentum 
that is relatively high, boosting the value of IZj with respect 
to that expected in a simple cooling model. 

Although plausible, we caution that the above discus- 
sion regarding the potential impact of the cold-mode to 
hot-mode accretion regime is highly speculative. In fact, 
the high resolution, hydrodynamical simulations of Sales 
et al. (2009), which in principle should naturally account 
for the cold mode vs. hot mode effects mentioned above, did 
not yield (disc) galaxies with the correct relation between 
£gf and IZj. Amazingly enough, their simulated galaxies 
revealed a tight, linear correlation between egf and IZj in 
perfect agreement with our predictions for a simple cool- 
ing model, despite the fact that their simulations included 
star formation and feedback (using several different imple- 
mentations). If taken at face value, these simulations seem 
to suggest that feedback has little to no effect on the egf 
vs. TZj relation; although it reduces 6gf, it causes a similar 
suppression of TZj. In addition, the Sales et al. simulations 
(which were only run to z = 2) seem to indicate that cold 
mode accretion, which should have been the dominant mode 
of accretion for the galaxies in their simulations, does not 
have the effect envisioned above. More detailed hydrody- 
namical simulations of galaxy formation (run all the way to 
z = 0) are needed to investigate these issues in more detail. 
In particular, we currently lack a proper understanding of 
how cold mode accretion and cold streams contribute to the 
build-up of mass and, in particular, angular momentum of 
disc galaxies. 
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